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We examine an exactly solvable model of decoherence - a spin-system interacting with a collec- 
tion of environment spins. We show that in this model (introduced some time ago to illustrate 
environment-induced superselection) generic assumptions about the coupling strengths lead to a 
universal (Gaussian) suppression of coherence between pointer states. We explore the regimes of 
validity of these results and discuss their relation to the spectral features of the environment and to 
the Loschmidt echo (or fidelity). Finally, we comment on the observation of such time dependence 
in spin echo experiments. 
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A single spin-system S (with states {|0),|1)}) inter- 
acting with an environment £ of many independent spins 
({Itfc} 7 II*)}) & = 1--ZV) through the Hamiltonian 

N 

H S£ = (|0) <0| - |1> <i|) f (Itfc) (1*1 - IU) (h\) 

k=l 

(1) 

may be the simplest solvable model of decoherence in 
spin systems. It was introduced some time ago Q to 
show that relatively straightforward assumptions about 
the dynamics can lead to the emergence of a preferred 
set of pointer states due to environment-induced super- 
selection (einselection) Q, [5] . Such models have gained 
additional importance in the past decade because of their 
relevance to quantum information processing 

The purpose of our paper is to show that - with a 
few additional natural and simple assumptions - one can 
evaluate the exact time dependence of the reduced den- 
sity matrix, and demonstrate that the off-diagonal com- 
ponents display a Gaussian (rather than exponential) de- 
cay. In effect, we exhibit a simple soluble example of 
a situation where the usual Markovian |j] assumptions 
about the evolution of a quantum open system are not 
satisfied at all times. Apart from their implications for 
decoherence, our results are also relevant to quantum er- 
ror correction |5( where precise precise knowledge of the 
dynamics is essential to select an efficient strategy. 

To demonstrate the Gaussian time dependence of de- 
coherence we first write down a general solution for the 
model given by Eq. Q). Starting with: 

N 

|*5fi(0)> = (a |0> + b |1» (a fe Itfc) + h || fe )) , (2) 

k=l 

the state of S£ at an arbitrary time is given by: 

|*5£(*)>=o|0>|£b(t)>+6|l)|f 1 (t)> (3) 

where 

N 

!$>(*)) - (g) (^ k e^ t/2 \h) + he~^ t/2 |4fe)) 
fc=l 

= M-t)). (4) 



The reduced density matrix of the system is then: 

Ps = Tre\*se(t)){9s£(t)\ 
= |a| 2 |0)(0|+a6V(t)|0)(l| 
+ a*6r*(t)|l)(0| + |fe| 2 |l)(l|, (5) 

where the decoherence factor r{t) = (£x(t)\£o(t)) can be 
readily obtained: 

N 

r(t) = XKWrfe^ + Wrfe-^ 1 ). (6) 

k=l 

It is straightforward to see that r(Q) = 1 and for t > 
it decays rapidly to zero, so that the typical fluctuations 
of the off-diagonal terms of ps will be small for large 
environments, since: 

(|r(t)| 2 )=2^n( 1 + (K| 2 -|^| 2 ) 2 )' ( 7 ) 
fc=l 

Here (...) denotes a long time average lj. Clearly, 
(l r W| 2 ) — '' 0, leaving p$ approximately diagonal in 

N^oc 

a mixture of the pointer states {|0},|1}} which retain 
preexisting classical correlations. 

This much was known since . The aim of this paper 
is to show that, for a fairly generic set of assumptions, 
the form of r(t) can be further evaluated and that - quite 
universally - it turns out to be approximately Gaussian 
in time. Thus, the simple model of Ref. Q predicts a uni- 
versal (Gaussian) form of the loss of quantum coherence, 
whenever the couplings gk of Eq. ^ are sufficiently con- 
centrated near their average value so that their standard 
deviation Ug^ — {gk)) 2 ) exists and is finite. When this 
condition is not fulfilled other sorts of time dependence 
become possible. In particular, r(t) may be exponential 
when the distribution of couplings is a Lorentzian. 

To obtain our main result we carry out the multiplica- 



tion of Eq.JfjJ), re-expressing r(t) as a sum: 



r(i) = J[\a k \ 2 e 
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J] |a fc | 2 e [ lt (- 9 '- 9 '" +E ^- 9 ")] + ... (8) 



There are 



etc. terms in the consecutive 



(o)> (7)> (2) 

sums above. The binomial pattern is clear, and can be 
made even more apparent by assuming that a k = a and 
gk = g for all k. Then, 
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r(t) 
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(9) 
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binomial 

N 



expansion of r(t) 



(\a\ 

We now note that, as follows from the Laplace-de 
Moivre theorem |(|, for sufficiently large TV the coeffi- 
cients of the binomial expansion of Eq. 10 can be ap- 
proximated by a Gaussian: 



2(N-l) l 



exp 



(1-N\f3\ 2 



2N\af3\- 



(10) 



This limiting form of the distribution of the eigenenergies 
of the composite S£ system immediately yields our main 
result: r(t) is approximately Gaussian since it is a Fourier 
transform of an approximately Gaussian distribution of 
the energies resulting from all the possible combinations 
of the couplings with the environment. 

The set of all the resulting energies must have an (ap- 
proximately) Gaussian distribution. This behavior is 
generic, a result of the law of large numbers 0: these 
energies can be thought of as being the terminal points 
of an TV-step random walk. The contribution of the A;-th 
spin of the environment to the random energy is +g or 
— g with probability \a\ 2 or |/3| 2 respectively (Fig. 1 a). 

The same argument can be carried out in the more 
general case of Eq. ©. The "random walk" picture that 
yielded the distribution of the couplings remains valid 
(see Fig. 1 b). However, now the individual steps in the 
random walk are not all equal. Rather, they are given 
by the set {gk} (see Eq. ^| with each step gk taken just 
once in a given walk. There are 2 N such distinct ran- 
dom walks, each contributing with the weight given by 
the product of the relevant \otk\ 2 and |/3fc| 2 to the sum 
of Eq. (JHJ). This exponential proliferation of the con- 
tributing coupling energies allows one to anticipate rapid 
convergence to the universal Gaussian form of the deco- 
herence factor r(t). 




-9k) 

3 k \ 2 as in Eq. JT^ 
2 = 1/2). (a) When all the couplings have the same size 
g (Eq. ©), a simple Newton's triangle leads to an ap- 



FIG. 1: The distribution of the energies obtains from the ran- 
dom walks with the steps given by the coupling size and in 
the direction (+gk or —gk) biased by the probabilities jo*] 2 
and \Pk\ 2 as in Eq. 115H (although in these examples we set 
\a k 
9k 

proximate Gaussian for the distribution of energies, (b) When 
the couplings differ from step to step (Eq. @), the resulting 
distribution still has a approximately Gaussian envelope for 
large N. 



Indeed, we can regard the energies resulting from the 
sums of gfc's as a random variable. Its probability distri- 
bution is given by products of the corresponding weights. 
That is, the typical term in Eq. [HJis of the form: 
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The resulting terminal energy is 



Ew = 



= E gk ~ E 5fe ' 

k£W+ kew- 



(11) 



(12) 



and the cumulative weight p\y is given by the correspond- 
ing product of \ak\ 2 and \Pk\ 2 - Each such specific ran- 
dom walk W corresponding to a given combination of 
right (k S W + ) and left (k <E W~) "steps" (see Fig. 1) 
contributes to the distribution of energies only once. The 
terminal points Ew may or may not be degenerate: As 
seen in Fig. 1, in the degenerate case, the whole collec- 
tion of 2 N random walks "collapses" into N + 1 terminal 
energies. More typically, in the degenerate case (also dis- 
played in Fig. 1), there are 2 N different terminal energies 
Ew In both cases, the "envelope" of the distribution 
P(Eyy) should be Gaussian, as we shall argue below. 

We note that the decoherence factor r(f) can be viewed 
as the characteristic function jj] (i.e., the Fourier trans- 
form) of the distribution of eigenenergies Ew- Thus, 



it) 



e iEt V (E)dE, 



(13) 
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FIG. 2: (Left panels) Assumed distribution of the couplings 
gt ■ (Center panels) Resulting distribution of the eigenenergies 
Ew (center panels) for TV = 6 (E w < 0) and TV = 24 (E w > 
0). In the case of |a*| 2 = 1/2 this distribution is in effect the 
"strength function" (local density of states). (Right panels) 
Decoherence factor r(t) for different initial conditions with 
TV = 6 (dashed lines), TV = 24 (thin solid lines) and the 
average (bold line). 

where the strength function i](E), also known as the local 
density of states (LDOS) Q is defined in general as 

r 1 (E) = '£\(*ss(0)\<Px)\ 2 S(E-E x ). (14) 

A 

Above \4>x) are the eigenstates of the full Hamiltonian 
and E\ its eigenenergies. In our particular model (Eq.QJ 
the eigenstates are associated with all possible random 
walks in the set W, and therefore 

r,(E) =J2pwS(E - E w ). (15) 
w 

The discussion of decoherence in our model is thus di- 
rectly related to the study of the characteristic function 
of the distribution of coupling energies rj(E). Moreover, 
since the E w 's are sums of gk's (that we assume inde- 
pendent of each other), r(t) is itself a product of char- 
acteristic functions of the distributions of the couplings 
{gk}, as we have already seen in the example of Eq. 
Thus, the distribution of Ew belongs to the class of the 
so-called infinitely divisible distributions 

The behavior of the decoherence factor r(t) - char- 
acteristic function of an infinitely divisible distribution 
- depends only on the average and variance of the dis- 
tributions of couplings weighted by the initial state of 
the environment |fj, |8| ■ The remaining task is to calcu- 
late rj(E), which can be obtained through the statistical 
analysis of the weighted random walk picture described 
above. If we denote Xk the random variable that takes 
the value or — gk with probability \ctk\ 2 or \(3k\ 2 re- 



spectively, then its mean value a& and its variance bk are 

ak = (\ak\ 2 - \0k\ 2 )gk, 

bl = g 2 k -al = A\ ak \ 2 \h\ 2 g 2 k. (16) 

The behavior of the sums of TV random variables Xk (and 
thus, of their characteristic function) depends on whether 
the so-called Lindeberg condition holds 6]. It is ex- 
pressed in terms of the cumulative variances B 2 ^ — ^ i>|, 
and it is satisfied when the probability of the large indi- 
vidual steps is small; e.g.: 

n ma* k -a k \> tB n ) — » 0, (17) 

l<k<N N^oo 

for any positive constant r. In effect, Lindeberg condition 
demands that B n be finite: when it is met, the resulting 
distribution of energies E = ^Xk is Gaussian 

P ( < — f (18) 
where En = ^ fc ak- In terms of the LDOS this implies 

^~-v^M =s ^ st )' (19) 

an expression in excellent agreement with numerical re- 
sults already for modest values of TV. This distribution of 
energies yields a corresponding approximately Gaussian 
time-dependence of r(t), as seen in Fig. 2. Moreover, 
at least for short times of interest for, say, quantum er- 
ror correction, r(t) is approximately Gaussian already 
for relatively small values of TV. This conclussion holds 
whenever the initial distribution of the couplings has a 
finite variance. The general form of r(t) after applying 
the Fourier transform of Eq. i|13|) is 

r(i) ~ e m »*e- B »*I 2 . (20) 

It is also interesting to investigate cases when Lin- 
deberg condition is not met. Here, the possible limit 
distributions are given by the stable (or Levy) laws 0. 
One interesting case that yields an exponential decay of 
the decoherence factor corresponds to the case of the 
Lorentzian distribution of couplings (see Fig. 3). Fur- 
ther intriguing questions concern the robustness of our 
conclusion under the changes of the model. We shall ad- 
dress this issue elsewhere 9] but, for the time being, we 
only note that the addition of a strong self-Hamiltonian 
proportional to a x changes the nature of the time decay 
jlfjj. On the other hand, small changes of the environ- 
ment Hamiltonians (like for instance dipolar interactions) 
seem to preserve the Gaussian nature of r(t). 

It is interesting to notice that the Fourier transform 
of the strength function n(E) is also related to the 
Loschmidt echo ^]] (or fidelity) in the so called Fermi 
Golden rule regime. The fact that the purity and the 
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FIG. 3: Same as Fig. 2 but for a Lorentzian distribution of 
the couplings gt- In this case r(t) decays exponentially. The 
histogram and the dashed line in r(t) correspond to N = 20, 
the straight thin line is a particular case for N = 100 and 
the thick line is the average. We note that the convergence 
is slower than in the Gaussian case of Fig. 2, because real- 
izations of gk are more likely to have one or two dominant 
couplings. Therefore, although the average shows a clear ex- 
ponential decay, fluctuations are noticeable even for large N. 
Notice also that the logarithmic scale confirms the long time 
saturation of r(t) at ~ 2 _JV/2 , Eq (jj. 



fidelity have closely related decay rates has been recently 
shown f° r the case of a bath composed of non- 
interacting harmonic oscillators. In this sense our results 
could be interpreted as an extension of the discussion of 
Ref. 12] to spin environments. 

The connection with fidelity is more easily seen if we 
write a generalized version of the Hamiltonian (|T|). 

«S£ = i(|0><0|®Wg + |l)<l|®?4). (21) 

The decoherence factor is then the overlap of the initial 
state of the environment |\E , g(0)) evolved with two differ- 
ent Hamiltonians, 

r(t) = <* £ (0)| e m °£ t/2 e- m i t/2 |*g(0)> , (22) 

which clearly has the form of the amplitude of the 
Loschmidt echo for the environment with the two states 
of the system as the perturbation. In the particular 
model that we are treating, 7i° £ = —H^ and thus 



r(t) = (* £ (0)|e- 



-mit 



*e(0)> 



(23) 



This expression is the survival probability of the initial 
state of the environment under the action of the Hamilto- 
nian Tig, which has been shown to be the Fourier trans- 
form of the strength function |l3| . This connection pro- 
vides another way to understand Eq. (IT51) 

Possible experimental applications of our considera- 
tions are in nuclear magnetic resonance, but also in other 
situations where two-level systems interact with spin en- 
vironments. We note that a Gaussian time dependence 
has been seen in the NMR setting [l4| but it is usually 
explained by spin diffusion models (which have rather 
different character and employ a different set of assump- 
tions). Moreover, there is a substantial body of work 
|lfjl Hil Hfi ] on decoherence due to spin environments, 



stimulated in part by the interests of quantum compu- 
tation. The relation between the decoherence factor and 
the strength function might prove useful in the physi- 
cal setting of strongly interacting fermions, where it has 
been shown that the strength function takes a Gaussian 
shape 01 ■ It is our hope that the simple analytic model 
described here will assist in gaining further insights into 
these fascinating problems. 
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